Single-cycle megawatt terahertz pulse generation from a wavelength-scale plasma 

oscillator driven by ultrashort laser pulses 
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The tremendous applications of terahertz (THz) spectroscopy and imaging require THz sources 
in different parameters. We propose a novel scheme to generate single-cycle powerful THz pulses by 
ultrashort intense laser pulses incident obliquely on a tenuous plasma slab of few THz wavelengths 
in thickness. This is made possible by driving a large amplitude electron plasma wave in the plasma 
slab, thus producing a net transient current at the plasma surfaces. Theory and simulations show 
that such a THz source is capable of providing megawatt power and field strength of MV/cm, which 
may open up new horizons for nonlinear THz science and applications. 

PACS numbers: 42.65.Re, 42.72.Ai, 52.25. Os, 52.59. Ye 
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Terahertz spectroscopy can probe the spectral prop- 
erties of molecules in a previously inaccessible electro- 
magnetic spectrum 1]. Its applications include the char- 
acterization of semiconductors 0] and high-temperature 
superconductors [|[, T-ray imaging of biomedical tis- 
sues d|, cellular structures @ and dielectric substances 
and the manipulation of bound atoms 0]. Most 
applications are based on the techniques of THz time- 
domain spectroscopy (THz-TDS) [if, where typically co- 
herent broadband THz pulses in the 2-5 THz frequency 
bandwidth are employed. Most broadband pulsed THz 
sources are based on the excitation of different materials 
with ultrashort laser pulses pj, such as through photo- 
conduction or optical rectification, but the output power 
is limited by the damage threshold of the optical materi- 
als used. 

So far, there have been no development in the non- 
linear regime due to the lack of sufficiently high power 
THz sources [lj. To address this problem a few efforts 
have been made using ultrashort electron bunches pro- 
duced by lasers or accelerators through the mechanisms 
of transition radiation and/or synchrotron radiation Q. 

Recently, by four- wave mixing in air with a laser inten- 
sity of ^10 14 W/cm 2 , single-cycle THz pulses with a field 
strength greater than kV/cm can be produced @, after 
the laser field-ionization and plasma generation thresh- 
old is surpassed. Whether this kind of THz emission will 
saturate for light intensities >10 15 W/cm 2 is an open 
question. 

For light intensities ^10 14 W/cm 2 , the leading part 
of the laser pulse can completely ionize the material and 
the interaction process becomes a pure laser-plasma in- 
teraction. Plasma has no thermal damage threshold and 
can sustain extremely intense light. The peak intensity 
of lasers today can be as high as 10 20 W/cm 2 and in the 
future is expected to reach 10 23 W/cm 2 Exploring 
new THz emission mechanisms in the context of intense 
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laser plasma interactions may produce higher power THz 
sources. It has been shown that laser wakefields (elec- 
tron plasma waves driven by the ponderomotive force 
of a laser pulse) in inhomogeneous plasmas can radiate 
high-efficiency THz waves at high power through linear 
mode conversion [ll|, [l^] • The THz pulses produced in 
this way are generally in multi-cycles and have a negative 
frequency chirp. 

The present work introduces a new THz emission 
mechanism in laser plasma interaction, which can di- 
rectly generate a single-cycle THz pulse with a field 
strength comparable with our earlier proposal [TUtHI- I n 
addition to the above applications, the single-cycle pulse 
has special implications for THz propagation physics and 
seismic surveys [Hj]. In the one-dimensional (ID) case, 
an electron plasma wave [Tij in a uniform plasma is de- 
scribed by Sn — 5n p exp[i(k p x — u> p t)], where 5n p is the 
density perturbation amplitude, lo p = y/ ne 2 /mso is the 
plasma frequency of the background plasma of density 
n, — e and m are the electron charge and mass respec- 
tively, and k p — lu p /c — 2t:/X p , where X p is the plasma 
wavelength. This infinite plasma wave can never emit 
electromagnetic waves at frequency u> p , since its dis- 
placement current (eodE/dt) exactly compensates the 
plasma current (—era;). However, for a finite plasma 
wave of length L ~ X p , there could be some electro- 
magnetic radiation. Firstly, for such a few-wavelength 
plasma oscillator, its displacement and plasma currents 
cannot completely counteract each other, in particular 
near the plasma boundaries. Secondly, since the plasma 
skin depth of the radiation at frequency u> p is k~ l , which 
is comparable to the plasma length L, the radiation can 
tunnel through the plasma into vacuum. Figure 1(a) is a 
schematic of this THz wave emission mechanism. 

In the following we give a theoretical analysis of this 
THz emission mechanism. For the interaction geome- 
try shown in Fig. 1(a), by Lorentz transformations we 
transform all the physical quantities from the laboratory 
frame to a moving frame of velocity csin^e^, where e y is 
the unit vector along the y direction. An electromag- 
netic wave with uj l = to', k L = (±fc' cos 8, k' sin 0, 0) 
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in the laboratory frame becomes oj m — uj' cos 9, k M = 
(±fc' cos 9, 0, 0) in the moving frame, wherein all the elec- 
tromagnetic waves propagate along the ±x directions, 
i.e. THz emission in the laboratory frame must be in 
the specular reflection and laser transmission directions. 
Plasma (electrons and ions) in the moving frame streams 
along — e y , with a relativistic factor 7 M = 1/ cos 6*. Fol- 
lowing Ref. [l5| and using the quasi-static approximation 
[l6|, we obtain the coupled equations (in SI units) 



dx 2 



„2 Q t 2 



-fs(x,i), 



5n = 



1 



d 2 <P co 2 
-7-17 = -irbn 

ox z & 



1 + (&l — e y tan( 



2cos6» 



(l/cos0 



- 1 



(1) 



(2) 



(3) 



where s(x,t) = —Sntan9e y / r f is the THz radiation 
source, &t and are the respective vector potentials 
of the THz wave and incident laser normalized by mc/e, 
4> is the scalar potential of the driven plasma wave nor- 
malized by mc 2 /e, Sn is the density perturbation of the 
plasma wave normalized by the initial plasma density n, 
7 = yjl + (sll — t&n9e y ) 2 + p 2 is the electron relativis- 
tic factor, and p x is the electron longitudinal momentum 
normalized by mc. 

The generation of the THz radiation is determined by 
Eq. (1), from which the electric field in the laboratory 
frame is found to be [151 
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where L is the plasma length, u> the laser frequency, 
and normalization with respect to mujc/e has been per- 
formed. Equation (4) shows that the THz emission is al- 
ways p-polarized. In the weakly relativistic case <j> « 1, 
Eqs. (2) and (3) lead to Sn oc a 2 L cos 9. In the THz emis- 
sion process, we assume 7 ~ 7 M . Substituting these into 
Eq. (4) we obtain 
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. p . From Eq. (5) we can see that there is no 



for L ~ A r 

THz emission for normal incidence (9 — 0), since there is 
no transverse electric field component. 

To test our proposal and the scaling rule in Eq. (5) 
we conduct a series of ID particle-in-cell (PIC) simula- 
tions. Taking into account the oblique incidence of the 
laser beam, our ID-PIC code adopts a moving frame as 
discussed above |15| and outputs all physics quantities 
in the laboratory frame. The initial plasma density is 
taken to be n = 0.0001n c , where n c = meoOj 2 /e 2 = 
1.1 x 10 21 (/xm/A) 2 cm~ 3 is the critical density for the 
laser pulse of wavelength A in vacuum. For A = 1 /xm, 



n = 1.1 x 10 17 cm -3 . The corresponding plasma fre- 
quency is uj p /2ir = 2.98 THz, which represents the central 
frequency of the THz emission. The plasma wavelength 
is A p = y/n c /nX = 100A. The incident laser pulse has a 
sine-square profile = eE^jmioc — ao sin 2 [7r(x— ct)/di] 
for < x — ct < dL, where is the laser pulse du- 
ration. Here ao is related to the peak laser intensity 
through 1 = a% - 1.37 x 10 18 (^m/A) 2 W/cm 2 . The rel- 
ativistic intensity threshold is reached at ao = 1. The 
laser pulse enters the left boundary of the simulation box 
with s polarization in order to distinguish it easily from 
the p-polarized THz emission from the wakefield. For the 
sine-square laser pulse, the excited wakefield amplitude is 
maximum when c?i w A p 12]. The simulation results also 
confirm that the THz emission is strongest for d^ = A p , 
thus we always set c?l = A p in the following. 

Defining fields F± = (E y ± cB z )/2 in the moving 
frame, we see that F + and F^. represent the forward 
and backward p-polarized electromagnetic waves, respec- 
tively. Tracing F + and F_ at the right and left bound- 
aries of the simulation box, we can obtain the tempo- 
ral profile of the radiated THz pulses in the reflection 
and transmission directions. Figure 1(b) shows the peak 
field strengths |-F±| max of the THz pulses as a function 
of the plasma length L. The laser pulse parameters are 
ao = 0.5, dL = 100A and 9 = 45°. The incident laser in- 
tensity is about 3.4 x 10 17 W/cm 2 (A = 1 am). We find 
that the plasma length L should be within [0.25A P , 2.7A p ] 
for intense THz pulses to be generated. When L > 3A P , 
the THz pulse amplitude decreases dramatically. 

We also find that the radiated THz pulse is single- 
cycle. Figures 1(c) and 1(d) show the time evolution of 
field components E x and B z in the x space for L = 100A. 
The field E x in Fig. 1(c) includes the longitudinal field 
of the wakefield and the electric field of the p-polarized 
THz emission. The wakefield is completely localized in 
the plasma region, while the electric field of the THz 
wave is mainly outside the plasma slab. In Fig. 1(d) B z 
is the pure magnetic field of the THz wave. It is obvious 
that two single-cycle THz pulses are radiated from the 
plasma region. Due to the propagation delay of the laser 
pulse, the pulse in the backward (reflection) direction is 
generated earlier than that in the forward (transmission) 
direction. For the specific laser wavelength A = 1 fj,m 
in Fig. 1(d), the field strength is found to be above 10 
MV/cm. When L > 3A P , the THz pulse is no longer 
single-cycle. 

Figures 2(a) and 2(b) illustrate the temporal profiles 
of the THz pulses shown in Figs. 1(c) and 2(d) together 
with two other incident angles of 30° and 60°. The shape 
of the transmitted THz wave F + is the same as that of 
the reflected F_. For 9 = 30°, the THz pulses have 
two cycles. Single-cycle THz emission is produced when 
9 > 45°. With increasing incident angle, the number of 
cycles included in the THz pulse decreases. Figure 2(c) 
displays the power spectra of the THz pulses F_. The 
central frequency is at 3 THz (A = 1 ^tm). The spec- 
trum width increases with the incident angle, because of 
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the shorter THz duration for the larger 9. The band- 10425416, 10390160, 10476033, and 10576035), the Na- 
widths approach 3-6 THz, meeting the requirements for tional High-Tech ICF Committee in China, and the 
the THz-TDS system. Figure 2(d) shows that the peak Knowledge Innovation Program, CAS. 
field strengths |F±| max are proportional to sin#, which 
agrees with Eq. (5). There is no THz emissions for 8 = 0. 

Figures 3(a) and 3(b) show that the THz field strength 
is proportional to both the laser intensity, i.e. a 2 , and the 
square root of the plasma density n. At lower intensities 
of 10 14 -10 15 W/cm 2 , the THz field strength is several 
tens of kV/cm, comparable to that generated through 
four- wave mixing in air Q under the same intensities. 

The above ID theoretical analysis and ID-PIC simu- 
lations are valid as long as the laser spot size is large 
compared with the plasma wavelength. For a laser beam 
with a Gaussian profile exp(—r 2 /w L ) in transverse space, 
the ID model applies for vol 3> X p . It is easily under- 
stood that, for wl < A p , the radiation source size is 
smaller than the radiated wavelength, so that the gen- 
erated THz wave will diffract dramatically. In order to 
have collimated THz emission, roj, > A p should be main- 
tained. 

To illustrate the multi-dimensional properties of the 
THz emission, Fig. 4 shows the THz emission in a 2D- 
PIC simulation. It is found that there are indeed two 
single-cycle THz pulses in the reflection and transmission 
directions. For A = 1 /im, the THz field strength is 42 
MV/cm, i.e. an intensity of 2.5 x 10 12 W/cm 2 . Since 
the radiation radius is about 30^m, this is equivalent to 
a peak power of 70 MW. 

A uniform plasma slab several tens of microns long can 
be readily formed from the gas jet targets commonly used 
in high- harmonic generation experiments [l7| . Mean- 
while, our numerical simulations show that THz pulses 
emitted from a nonuniform plasma with a trapezoid den- 
sity profile are similar to those from a uniform plasma 
slab, provided that the ascending and decending parts of 
the trapezoid are also at the THz wavelength scale. This 
suggests that the density homogeneity of the plasma slab 
is not necessary. 

To conclude, we have presented a method for pro- 
ducing single-cycle high power THz radiation from 
wavelength-scale plasma oscillators (L ~ A p ). It is emit- 
ted by the transient net currents induced at the plasma 
surfaces while building-up the plasma oscillators. This 
mechanism together with that for THz emission by lin- 
ear mode conversion in inhomogeneous plasmas (L 3> A p ) 
[111 [l2j provide a complete picture for interpreting the 
early experimental observation of THz emission in in- 
tense laser plasma interaction |18j . 
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Figure Captions 

FIG. 1 (color online), (a) Schematic of THz emission 
from a few-wavelength long plasma oscillator excited by 
a laser pulse obliquely incident on a uniform plasma slab. 
The plasma length is L ~ \ p and 9 is the laser incident 
angle. THz waves are radiated in the specular reflec- 
tion and laser transmission directions, similar to the op- 
tical rectification scheme, (b) The peak field strengths 
|-F±|max of the radiated THz pulses as a function of the 
plasma length L. The plasma is of density n = 0.0001n c . 
The laser pulse parameters are ao = 0.5, = 100A and 



9 = 45°. Spatial-temporal plots of (c) the electric field E x 
and (d) the pure THz magnetic field B z for the plasma 
length L = 100A (located in 100A < x < 200A). The ver- 
tical dashed lines represent the plasma boundaries. The 
dashed arrow line marks the laser propagation trajectory. 
Other parameters are the same as in (b). 

FIG. 2 (color online). Temporal profile, frequency 
spectrum and field amplitude of the THz emission, (a) 
Temporal profiles of the THz waves F + for incident an- 
gles of 9 = 30°, 45° and 60°. (b) Temporal profiles of 
the THz waves Other parameters are the same as 
in Fig. l(c,d). (c) The power spectra of F_ in (b). (d) 
The peak field strengths |-F±| m ax of the THz pulses as a 
function of sin(9 with 9 G [0°, 15°, 30°, 45°, 60°, 75°]. 

FIG. 3 (color online). The scaling rule of the THz 
emission, (a) The peak field strengths \F± | max of the THz 
pulses as a function of the laser intensity. The plasma 
parameters are n = 0.0001n c , L = 50A and the laser 
pulse parameters dh = 100A and 9 = 45°. (b) |F±| max 
of the THz pulses as a function of the electron density 
n. The laser pulse parameters are ao = 0.5 and 9 = 45°. 
For a given plasma density n, we take L = di = A p . The 
dashed lines are fitted curves. 

FIG. 4 (color online). 2D spatial plot of the pure 
THz magnetic field B z from 2D-PIC simulation. We 
take n = 0.0025n c , corresponding to A p = 20A and 
u) p /2tt = 14.9 THz. The plasma length is L = 25A. The 
laser pulse is s-polarized, focused on the plasma slab sur- 
face, and has parameters a = 0.5, wl = 30A, d^ = 20A 
and 9 = 50°. The dashed rectangle shows the plasma 
region, the solid arrow marks the laser propagation axis, 
and the dashed arrows the THz emission directions, one 
along the laser propagation and another along the spec- 
ular reflection direction. 
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